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Abstract. One of the major advances needed to realize all-optical information 
processing of light is the ability to delay or coherently store and retrieve optical 
information in a rapidly tunable manner. In the classical domain, this optical 
buffering is expected to be a key ingredient to managing the flow of information over 
complex optical networks. Such a system also has profound implications for quantum 
information processing, serving as a long-term memory that can store the full quantum 
information contained in an optical pulse. Here we suggest a novel approach to light 
storage involving an optical waveguide coupled to an optomechanical crystal array, 
where light in the waveguide can be dynamically and coherently transferred into long- 
lived mechanical vibrations of the array. Under realistic conditions, this system is 
capable of achieving large bandwidths and storage/delay times in a compact, on-chip 
platform. 
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1. Introduction 

Light is a natural candidate to transmit information across large networks due to its 
high speed and low propagation losses. A major obstacle to building more advanced 
optical networks is the lack of an ail-optically controlled device that can robustly delay 
or store optical wave-packets over a tunable amount of time. In the classical domain, 
such a device would enable all-optical buffering and switching, bypassing the need to 
convert an optical pulse to an electronic signal. In the quantum realm, such a device 
could serve as a memory to store the full quantum information contained in a light pulse 
until it can be passed to a processing node at some later time. 

A number of schemes to coherently delay and store optical information are 
being actively explored. These range from tunable coupled resonator optical 
waveguide (CROW) structures [1, 2], where the propagation of light is dynamically 
altered by modulating the refractive index of the system, to electromagnetically induced 
transparency (EIT) in atomic media [3, 4], where the optical pulse is reversibly mapped 
into internal atomic degrees of freedom. While these schemes have been demonstrated 
in a number of remarkable experiments [5, 6, 7, 8], they remain difficult to implement 
in a practical setting. Here, we present a novel approach to store or stop an optical 
pulse propagating through a waveguide, wherein coupling between the waveguide and 
a nearby nano-mechanical resonator array enables one to map the optical field into 
long-lived mechanical excitations. This process is completely quantum coherent and 
allows the delay and release of pulses to be rapidly and ail-optically tuned. Our scheme 
combines many of the best attributes of previously proposed approaches, in that it 
simultaneously allows for large bandwidths of operation, on-chip integration, relatively 
long delay/storage times, and ease of external control. Beyond light storage, this work 
opens up the intriguing possibility of a platform for quantum or classical all-optical 
information processing using mechanical systems. 

2. Description of system: an optomechanical crystal array 

An optomechanical crystal [9] is a periodic structure that constitutes both a photonic [10] 
and a phononic [11] crystal. The ability to engineer optical and mechanical properties in 
the same structure should enable unprecedented control over light-matter interactions. 
Planar two-dimensional (2D) photonic crystals, formed from patterned thin dielectric 
films on the surface of a microchip, have been succesfuUy employed as nanoscale 
optical circuits capable of efficiently routing, diffracting, and trapping light. Fabrication 
techniques for such 2D photonic crystals have matured significantly over the last decade, 
with experiments on a Si chip [12] demonstrating excellent optical transmission through 
long (N > 100) linear arrays of coupled photonic crystal cavities. In a similar Si chip 
platform it has recently been shown that suitably designed photonic crystal cavities also 
contain localized acoustic resonances which are strongly coupled to the optical field via 
radiation pressure [9]. These planar optomechanical crystals (OMCs) are thus a natural 
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Figure 1. (a) Illustration of a double optical cavity system forming the unit cell 
of the optomechanical array. A two-way optical waveguide is coupled to a pair of 
optical cavity modes ai and a2, whose resonance frequencies differ by the frequency 
of the mechanical mode b. Both optical modes leak energy into the waveguide 
at a rate /^ex and have an inherent decay rate i^in- The mechanical resonator 
optomechanically couples the two optical resonances with a cross coupling rate of h. 
(b) A simplified system diagram where the classically driven cavity mode a2 is effectively 
eliminated to yield an optomechanical driving amplitude Qrn between the mechanical 
mode and cavity mode ai. (c) Frequency-dependent refiectance (black curve) and 
transmittance (red) of a single array element, in the case of no optomechanical driving 
amplitude Qrn = (dotted line) and an amplitude of Clrn = (solid line). The 

inherent cavity decay is chosen to be A^in = O.lA^ex- (inset) The optomechanical coupling 
creates a transparency window of width '^4r^^//^ex for a single element and enables 
perfect transmission on resonance, = 0. (d) Energy level structure of simplified 
system. The number of photons and phonons are denoted by ni and n^, respectively. 
The optomechanical driving amplitude Qrn couples states -h l^m^^i + 1) 

while the light in the waveguide couples states |n^,ni) l^m^^i + !)• The two 
couplings create a set of A-type transitions analogous to that in EIT. 

In the following we consider an optomechanical crystal containing a periodic array of 
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such defect cavities (see figures l(a),(b)). Each element of the array contains two optical 
cavity modes (denoted 1,2) and a co-localized mechanical resonance. The Hamiltonian 
describing the dynamics of a single element is of the form 

^om = huia\ai + hjuj2Ci\a2 + hcUrnb^b -\- hh (h -\- b^^ (^d\a2 + aldi^ . (1) 

Here a;i^2 are the resonance frequencies of the two optical modes, Um is the mechanical 
resonance frequency, and ai,a2,& are annihilation operators for these modes. The 
optomechanical interaction cross-couples the cavity modes 1 and 2 with a strength 
characterized by h and that depends linearly on the mechanical displacement xcc{b-\-b'^). 
While we formally treat ai,a2,6 as quantum mechanical operators, for the most part 
it also suffices to treat these terms as dimensionless classical quantities describing the 
positive-frequency components of the optical fields and mechanical position. In addition 
to the optomechanical interaction described by equation (1), the cavity modes 1 are 
coupled to a common two-way waveguide (described below). Each element is decoupled 
from the others except through the waveguide. 

The design considerations necessary to achieve such a system are discussed in detail 
in the section "Optomechanical crystal design." For now, we take as typical parameters 
uji/27r = 200 THz, Um/^n = 10 GHz, h/27r = 0.35 MHz, and mechanical and (unloaded) 
optical quality factors of = Cl;^/7^^10^ (room temperature)-10^ (low temperature) 
and Qi = uJi/tvi^in = 3x10^, where 7^ is the mechanical decay rate and tvi^in is 
the intrinsic optical cavity decay rate. Similar parameters have been experimentally 
observed in other OMC systems [9, 13]. In practice, one can also over-couple cavity 
mode 1 to the waveguide, with a waveguide-induced optical decay rate that is much 
larger than /^in. 

For the purpose of slowing light, the cavity modes 2 will be resonantly driven by 
an external laser, so that to good approximation a2^Q^2(^)e~*^^* can be replaced by its 
mean-field value. We furthermore consider the case where the frequencies are tuned such 
that (jJi = (jj2 -\-(jOm- Keeping only the resonant terms in the optomechanical interaction, 
we arrive at a simplified Hamiltonian for a single array element (see figure 1(b)), 

i/om = hujid\di + humb'^b + hQrnit) (at 6e"^^"^^""^-)^ + /i.e.) . (2) 

Here we have defined an effective optomechanical driving amplitude ^mit) = ha2{t) 
and assume that 0^2 (^) is real. Mode 2 thus serves as a "tuning" cavity that mediates 
population transfer (Rabi oscillations) between the "active" cavity mode 1 and the 
mechanical resonator at a controllable rate i^(t), which is the key mechanism for our 
stopped-light protocol. In the following analysis, we will focus exclusively on the active 
cavity mode and drop the "1" subscript. 

A Hamiltonian of the form (2) also describes an optomechanical system with a 
single optical mode, when the cavity is driven off resonance at frequency uji — ujjn and 
ai corresponds to the sidebands generated at frequencies ±uJm around the classical 
driving field. For a single system, this Hamiltonian leads to efficient optical cooling 
of the mechanical motion [14, 15], a technique being used to cool nano-mechanical 
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systems toward their quantum ground states [16, 17, 18, 19]. While the majority 
of such work focuses on how optical fields affect the mechanical dynamics, here we 
show that the optomechanical interaction strongly modifies optical field propagation 
to yield the slow/stopped light phenomenon. Equation (2) is quite general and thus 
this phenomenon could in principle be observed in any array of optomechanical systems 
coupled to a waveguide. In practice, there are several considerations that make the 2D 
OMC "ideal." First, our system exhibits an extremely large optomechanical coupling h 
and contains a second optical tuning cavity that can be driven resonantly, which enables 
large driving amplitudes Qrn using reasonable input power [20]. Using two different 
cavities also potentially allows for greater versatility and addressability of our system. 
For instance, in our proposed design the photons in cavity 1 are spatially filtered from 
those in cavity 2 [20]. Second, the 2D OMC is an easily scalable and compact platform. 
Finally, as described below, the high mechanical frequency of our device compared to 
typical optomechanical systems allows for a good balance between long storage times 
and suppression of noise processes. 



3. Slowing and stopping light 

3.1. Static regime 

We first analyze propagation in the waveguide when flmit) = is static during the 
transit interval of the signal pulse. As shown in the Appendix, the evolution equations 
in a rotating frame for a single element located at position zj along the waveguide are 
given by 



f = -^b + inm& + F^{t). (4) 

Equation (3) is a standard input relation characterizing the coupling of right- (ciR^i^) and 
left-propagating {aL,in) optical input fields in the waveguide with the cavity mode. Here 
K = /^ex + ^in is the total optical cavity decay rate, a^{z) is quantum noise associated 
with the inherent optical cavity loss, and for simplicity we have assumed a linear 
dispersion relation Uk = c\k\ in the waveguide. Equation (4) describes the optically 
driven mechanical motion, which decays at a rate 7^ and is subject to thermal noise 
Fj^it). The cavity mode couples to the right-propagating field through the equation 

(c^ + ^ni^^t) = ~ ^i)^ + ikoO^R, (5) 

where ko = cji/c. We solve the above equations to find the reflection and transmission 
coefficients r, t of a single element for a right-propagating incoming field of frequency 
Uk (see Appendix). In the limit where 7^ = 0, and defining 5k=uJk — cji, 
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while t = l-\-r. Example reflectance and transmittance curves are plotted in figure 1(c). 
For any non-zero fi^, a single element is perfectly transmitting on resonance, whereas for 
Qrn = resonant transmission past the cavity is blocked. When J^^t^O, excitation of the 
cavity mode is inhibited through destructive interference between the incoming field and 
the optomechanical coupling. In EIT, a similar effect occurs via interference between 
two electronic transitions. This analogy is further elucidated by considering the level 
structure of our optomechanical system (figure 1(d)), where the interference pathways 
and the "A" -type transition reminiscent of EIT are clearly visible. The interference 
is accompanied by a steep phase variation in the transmitted field around resonance, 
which can result in a slow group velocity. These steep features and their similarity to EIT 
in a single optomechanical system have been theoretically [21, 22] and experimentally 
studied [23, 24], while interference effects between a single cavity mode and two 
mechanical modes have also been observed [25]. From r, t for a single element, the 
propagation characteristics through an infinite array (figure 2(a)) can be readily obtained 
via band structure calculations [10]. To maximize the propagation bandwidth of the 
system, we choose the spacing d between elements such that k^d = (2n + l)7r/2 where 
n is a non-negative integer. With this choice of phasing the reflections from multiple 
elements destructively interfere under optomechanical driving. Typical band structures 
are illustrated in flgures 2(b)-(f). The color coding of the dispersion curves (red 
for waveguide, green for optical cavity, blue for mechanical resonance) indicates the 
distribution of energy or fractional occupation in the various degrees of freedom of the 
system in steady-state. Far away from the cavity resonance, the dispersion relation is 
nearly linear and simply reflects the character of the input optical waveguide, while the 
propagation is strongly modified near resonance (uu = oui = 002 -\- cjm)- In the absence 
of optomechanical coupling (Jl^ = 0), a transmission band gap of width ^ k forms 
around the optical cavity resonance (refiections from the bare optical cavity elements 
constructively interfere). In the presence of optomechanical driving, the band gap splits 
in two (blue shaded regions) and a new propagation band centered around the cavity 
resonance appears in the middle of the band gap. For weak driving (f]^</^) the width of 
this band is ^4f^^//^, while for strong driving (f^^>A^) one recovers the "normal mode 
splitting" of width ^2J~i^ [26]. This relatively fiat polaritonic band yields the slow-light 
propagation of interest. Indeed, for small Qrn the steady-state energy in this band is 
almost completely mechanical in character, indicating the strong mixing and conversion 
of energy in the waveguide to mechanical excitations along the array. 

It can be shown that the Bloch wavevector near resonance is given by (see 
Appendix) 



The group velocity on resonance, Vg = {dkQff/d6k)~^\sk=o = 2(ii7^//^ex can be 
dramatically slowed by an amount that is tunable through the optomechanical coupling 
strength J"^^. The quadratic and cubic terms in fceff characterize pulse absorption and 
group velocity dispersion, respectively. In the relevant regime where /^ex»^in, ^ex^^m. 



l^eyi^k ^^ex^in^^ (^^ex S/^ex^hi -'-^^ex^m)^^ 



(7) 
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Figure 2. (a) Illustration of an optomechanical crystal array. A two-way optical 
waveguide is coupled to a periodic array of optomechanical elements spaced by a 
distance d. The optical cavity modes aj of each element leak energy into the waveguide 
at a rate k^q^ and have an inherent decay rate Hii^. The mechanical resonator of each 
element has frequency and is optomechanically coupled to the cavity mode through 
a tuning cavity (shown in figure 1) with strength (b) The band structure of the 
system, for a range of driving strengths between = and = ^- The blue 
shaded regions indicate band gaps, while the color of the bands elucidates the fractional 
occupation (red for energy in the optical waveguide, green for the optical cavity, and 
blue for mechanical excitations). The dynamic compression of the bandwidth is clearly 
visible as ftrn 0. (c) Band structure for the case Qrn = /^/lO is shown in greater 
detail, (d) The fractional occupation for each band in (c) is plotted separately. It can 
be seen that the polaritonic slow-light band is mostly mechanical in nature, with a 
small mixing with the waveguide modes and negligible mixing with the optical cavity 
mode. Zoom-ins of figures (c) and (d) are shown in (e) and (f). 
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these effects are negligible within a bandwidth Acj^min 



Ari/3 



)). The 



second term is the bandwidth over which certain frequency components of the pulse 
acquire a 7r-phase shift relative to others, leading to pulse distortion. This yields a 
bandwidth-delay product of 



for static Q^n and negligible mechanical losses. When intrinsic optical cavity losses are 
negligible, and if one is not concerned with pulse distortion, light can propagate over 
the full bandwidth ^4Q^//^ of the slow-light polariton band and the bandwidth-delay 
product increases to AcjTdeiay^^ (see Appendix). On the other hand, we note that if we 
had operated in a regime where kod = 7rn, constructive interference in reflection would 
limit the bandwidth-delay product to AcjTdeiay^l, independent of system size. 

In the static regime, the bandwidth-delay product obtained here is analogous 
to CROW systems [2]. In the case of EIT, a static bandwidth-delay product of 
AcjTdeiay^VOD results, where OD is the optical depth of the atomic medium. This 
product is limited by photon absorption and re-scattering into other directions, 
and is analogous to our result AcjTdeiay ^ V ^^ex/ ^in in the case of large intrinsic 
cavity linewidth. On the other hand, when /^in is negligible, photons are never lost 
and reflections can be suppressed by interference. This yields an improved scaling 
AcjTdeiay^^^^^ or depending on whether one is concerned with group velocity 

dispersion. In atomic media, the weak atom-photon coupling makes achieving OD > 100 
very challenging [27]. In contrast, in our system as few as A^^IO elements would be 
equivalently dense. 

3.2. Storage of optical pulse 

We now show that the group velocity Vg{t) = 2dQ'^{t) / k.q^ can in fact be adiabatically 
changed once a pulse is completely localized inside the system, leading to distortion-less 
propagation at a dynamically tunable speed. In particular, by tuning Vg{t)^0^ the pulse 
can be completely stopped and stored. 

This phenomenon can be understood in terms of the static band structure of the 
system (figure 2) and a "dynamic compression" of the pulse bandwidth. The same 
physics applies for CROW structures [1, 28], and the argument is re-summarized here. 
First, under constant fi^, an optical pulse within the bandwidth of the polariton band 
completely enters the medium. Once the pulse is inside, we consider the effect of a 
gradual reduction in Qrnit)- Decomposing the pulse into Bloch wavevector components, 
it is clear that each Bloch wavevector is conserved under arbitrary changes of as it 
is fixed by the system periodicity. Furthermore, transitions to other bands are negligible 
provided that the energy levels are varied adiabatically compared to the size of the gap, 
which translates into an adiabatic condition |(rf/G?t)(Q^//^)|</^^. Then, conservation of 
the Bloch wavevector implies that the bandwidth of the pulse is dynamically compressed, 
and the reduction in slope of the polariton band (figure 2) causes the pulse to propagate 




(8) 
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at an instantaneous group velocity Vg{t) without any distortion. In the hmit that Q^^O^ 
the polaritonic band becomes flat and completely mechanical in character, indicating 
that the pulse has been reversibly and coherently mapped onto stationary mechanical 
excitations within the array. We note that since is itself set by the tuning cavities, its 
rate of change cannot exceed the optical linewidth and thus the adiabaticity condition 
is always satisfied in the weak-driving regime. 

The maximum storage time is set by the mechanical decay rate, ^1/7^. For realistic 
system parameters a;m/27r = 10 GHz and Qm = 10^, this yields a storage time of ^10 /iS. 
In CROW structures, light is stored as circulating fields in optical nano-cavities, where 
state of the art quality factors of Q^IO^ limit the storage time to ^1 ns. The key 
feature of our system is that we effectively "down-convert" the high-frequency optical 
fields to low- frequency mechanical excitations, which naturally decay over much longer 
time scales. While storage times of ^10 ms are possible using atomic media [29], their 
bandwidths so far have been limited to < 1 MHz [28]. In our system, bandwidths of 
^1 GHz are possible for realistic circulating powers in the tuning cavities. 



3.3. Imperfections in storage 

The major source of error in our device will be mechanical noise, which through the 
optomechanical coupling can be mapped into noise power in the optical waveguide 
output. In our system, mechanical noise emerges via thermal fiuctuations and 
Stokes scattering (corresponding to the counter-rotating terms in the optomechanical 
interaction that we omitted from Equation (2)). To analyze these effects, it suffices 
to consider the case of static fi^, and given the linearity of the system, no waveguide 
input (such that the output will be purely noise). For a single array element, the 
optomechanical driving results in optical cooling of the mechanical motion [14, 15], 
with the mechanical energy Em evolving as (see Appendix) 

dE 

= -Im {Em - f^mnth) " Fopt^^m + ^opt _^ IQ^2 + * 

The first term on the right describes equilibration with the thermal surroundings, 
where nth = (e^^^^/^^^^ — 1)"-^ is the Bose occupation number at the mechanical 
frequency and T5 is the bath temperature. The second (third) term corresponds to 
cooling (heating) through anti-Stokes (Stokes) scattering, with a rate proportional to 
Topt = ^^m/i^' The Stokes process is suppressed relative to the anti-Stokes in the 
limit of good sideband resolution /^/a;^«l. For an array of elements, a simple 
upper bound for the output noise power at one end of the waveguide is given by 
^^noise = (V^) (ropt^ss) ^(cji/cl;^) (/^ex//^) , whcrc Egg is the steady-state solution of 
Equation (9). The factor of 1/2 accounts for the optical noise exiting equally from 
both output directions, Fopt^'ss is the optically-induced mechanical energy dissipation 
rate, and /^ex/^ describes the waveguide coupling efficiency. The term uji/ujm represents 
the transduction of mechanical to optical energy and is essentially the price that one 
pays for down-converting optical excitations to mechanical to yield longer storage times 
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- in turn, any mechanical noise gets "up-converted" to optical energy (whereas the 
probability of having a thermal optical photon is negligible). In the relevant regime 
where ropt»7m, 



This noise analysis is valid only in the weak-driving regime (J7^</^) [14, 15]. The strong 
driving regime, where the mechanical motion acquires non-thermal character [32] and 
can become entangled with the optical fields [33], will be treated in future work. 

At room temperature, nth^fc^T^/ZiCL;^ is large and thermal noise will dominate, 
yielding a noise power of ^0.4 nW per element for previously given system parameters 
and A^ex/^^1- This is independent of Fopt provided that ropt»7m, which refiects the 
fact that all of the thermal heating is removed through the optical channel. For high 
temperatures, the thermal noise scales inversely with o;^, and the use of high-frequency 
mechanical oscillators ensures that the noise remains easily tolerable even at room 
temperature. 

Thermal noise in the high-frequency oscillator can essentially be eliminated 
in cryogenic environments, which then enables faithful storage of single photons. 
Intuitively, a single-photon pulse can be stored for a period only as long as the 
mechanical decay time ^ 7"^, and as long as a noise-induced mechanical excitation is 
unlikely to be generated over a region covering the pulse length and over the transit time 
"^deiay The latter condition is equivalent to the statement that the power P^i^^huJiAuj in 
the single-photon pulse exceeds Pnoise- While we have focused on the static regime thus 
far, when thermal heating is negligible, realizing Pph/^^noise^l in the static case in fact 
ensures that the inequality holds even when Tomtit) is time- varying. Physically, the rate 
of Stokes scattering scales linearly with Fopt while the group velocity scales inversely, 
and thus the probability of a noise excitation being added on top of the single-photon 
pulse is fixed over a given transit length. 

In a realistic setting, the optomechanical driving amplitude itself will 
be coupled to the bath temperature, as absorption of the pump photons in the 
tuning cavities leads to material heating. To understand the limitations as a 
quantum memory, we have numerically optimized the static bandwidth-delay product 
AcjTdeiay for ^ train of single-photon pulses, subject to the constraints Auu < 



min (2V2J^y VA^^ex/^in, 2{67ryWj{t^,^N^/^)) , Pph/Pnoise > 1, and 7mrdeiay < 1. As 



a realistic model for the bath temperature, we take T5 = Tq + x^l = -\- x{^m/hY-) 
where To is the base temperature and x^2 /iK is a temperature coefficient that 
describes heating due to pump absorption (see Appendix). Using Tq = 100 mK and 
Qm = 10^, we find (AcL;rdeiay)max^llO, which is achieved for parameter values A^^275, 




(10) 



Kex/27r~l.l GHz, and O™/27r~130 MHz. 
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Figure 3. (a) Top view of the proposed optomechanical crystal array, with the 
superlattice unit-cell of length d highlighted in the center. The unit cell contains two 
coupled L2 defect cavities (shaded in grey) with two side-coupled linear defect optical 
waveguides (shaded in red). The different envelope functions pertain to the odd and 
even optical cavity supermodes (green solid and dashed lines, respectively), the odd 
mechanical supermode (blue solid line), and the odd and even optical waveguide modes 
(red solid and dashed lines). The displacement field amplitude |Q(r)| of the mechanical 
mode and in-plane electric field amplitude |E(r)| of the optical mode are shown in (b) 
and (c), respectively, for a single L2 defect cavity, (d) Bandstructure of the linear- 
defect waveguide (grey) and the zone folded superlattice of the entire coupled-resonator 
system (red). The cavity mode (green) crosses the superlattice band at mid-zone, and 
the waveguide-cavity interaction is shown in more detail in the insets (e) and (a) for 
the even {ay = -hi) and odd {(jy = — 1) supermodes, respectively. 



4. Optomechanical crystal design 

A schematic showing a few periods of our proposed 2D OMC slow-light structure is given 
in figure 3. The structure is built around a "snowflake" crystal pattern of etched holes 
into a Sihcon slab [20]. This pattern, when implemented with a physical lattice constant 
of a = 400 nm, snowflake radius r = 168 nm, and snowflake width w = 60 nm (see 
figure 3(a)), provides a simultaneous phononic bandgap from 8.6 to 12.6 GHz and a 
photonic pseudo-bandgap from 180 to 230 THz (see Appendix). Owing to its unique 
bandgap properties, the snowfiake patterning can be used to form waveguides and 
resonant cavities for both acoustic and optical waves simply by removing regions of the 
pattern. For instance, a single point defect, formed by removing two adjacent holes (a 
so-called "L2" defect), yields the co-localized phononic and photonic resonances shown 
in figures 3(b) and (c), respectively. The radiation pressure, or optomechanical coupling 
between the two resonances can be quantified by a coupling rate, which corresponds 
to the frequency shift in the optical resonance line introduced by a single phonon of the 
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mechanical resonance. Numerical finite-element-method (FEM) simulations of the L2 
defect indicate the mechanical resonance occurs at Cc;^/27r = 11.2 GHz, with a coupling 
rate of 5'/27r = 489 kHz to the optical mode at frequency uJo/27r = 199 THz (free-space 
optical wavelength of Aq ^ 1500 nm). 

In order to form the double-cavity system described in the slow-light scheme above, 
a pair of L2 cavities are placed in the near- field of each other as shown in the dashed box 
region of figure 3(a). Modes of the two degenerate L2 cavities mix, forming supermodes 
of the double-cavity system which are split in frequency. The frequency splitting between 
modes can be tuned via the number of snowfiake periods between the cavities. As 
described in more detail in the Appendix, it is the optomechanical cross-coupling of the 
odd (E_) and even (E+) optical supermodes mediated by the motion of the odd parity 
mechanical supermode (Q_) of the double-cavity that drives the slow-light behaviour 
of the system. Since Q_ is a displacement field that is antisymmetric about the two 
cavities, there is no optomechanical self-coupling between the optical supermodes and 
this mechanical mode. On the other hand, the cross-coupling between the two different 
parity optical supermodes is large and given hy h = g/V2 = 27r(346 kHz). By letting 
ai, a2 and b be the annihilation operators for the modes E_, E+ and Q_, we obtain the 
system Hamiltonian of equation (1). 

The different spatial symmetries of the optical cavity supermodes also allow them 
to be addressed independently. To achieve this we create a pair of linear defects in the 
snowfiake lattice as shown in figure 3(a), each acting as a single-mode optical waveguide 
at the desired frequency of roughly 200 THz (see figure 3(d). Sending light down both 
waveguides, with the individual waveguide modes either in or out of phase with each 
other, will then excite the even or odd supermode of the double cavity, respectively. 
The waveguide width and proximity to the L2 cavities can be used to tune the cavity 
loading (see Appendix), which for the structure in figure 3(a) results in the desired 
/^ex/27r = 2.4 GHz. It should be noted that these hne-defect waveguides do not guide 
phonons at the frequency of Q_, and thus no additional phonon leakage is induced in 
the localized mechanical resonance. 

The full slow-light waveguide system consists of a periodic array of the double- 
cavity, double- waveguide structure. The numerically computed band diagram, for 
spacing d = 15a periods of the snowfiake lattice between cavity elements (the 
superlattice period), is shown in figure 3(d). This choice of superlattice period results 
in the folded superlattice band intersecting the E_ (ai) cavity frequency ui at roughly 
mid-zone, corresponding to the desired inter-cavity phase shift of kd = 7r/2. A zoom-in 
of the bandstructure near the optical cavity resonances is shown in figures 3(e) and (f). 
In figure 3(e) the even parity supermode bandstructure is plotted (i.e., assuming the 
even supermode of the double- waveguide is excited), whereas in figure 3(f) it is the odd 
parity supermode bandstructure. 

A subtlety in the optical pumping of the periodically arrayed waveguide system 
is that for the E+ (^2) optical cavity resonance at uj2^ there exists a transmission 
bandgap. In order to populate cavity a2, then, and to create the polaritonic band 
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Figure 4. The input power (red line) required to achieve the system parameters used 
in the text, i.e. l]^/27r = 130 MHz with h/27T = 0.346 MHz, and the attenuation per 
unit ceh a (sohd black line) are shown as a function of detuning of the pump beam 
from the pump cavity frequency. The dotted line is the approximate expression derived 
for the attenuation, a ^ ^ex^in/4(^^. The grey region indicates the band gap in which 
the pump cavities cannot be excited from the waveguide. The trade-off between small 
pump input powers and low pump attenuation factors is readily apparent in this plot. 

at cji, the pump beam must be slightly off-resonant from 6^2, but still at uji — uOm- 
We achieve this by choosing a double-cavity separation (14 periods) resulting in a 
cavity mode splitting ((c^i — = 9.7 GHz) slightly smaller than the mechanical 

frequency {um/"^^ = 11-2 GHz), as shown in figures 3(e) and (f). By changing the 
detuning between uji — uj2 and cj^, a trade-off can be made between the attenuation 
of the pump beam per unit cell, a = exp(— Im{i^}rf), and total required input power, 
shown in figure 4. In Appendix D.3 we show that the total attenuation per unit cell, is 
given by a ^ Interestingly, by using higher input powers such that a ^ 0, 

it is in principle possible to eliminate completely effects due to absorption in the array, 
which may lead to inhomogeneous pump photon occupations. 

5. Outlook 

The possibility of using optomechanical systems to facilitate major tasks in classical 
optical networks has been suggested in several recent proposals [35, 25, 20]. This present 
work not only extends these prospects, but proposes a fundamentally new direction 
where optomechanical systems can be used to control and manipulate light at a quantum 
mechanical level. Such efforts would closely mirror the many proposals to perform 
similar tasks using EIT and atomic ensembles [4]. At the same time, the optomechanical 
array has a number of novel features compared to atoms, in that each element can 
be deterministically positioned, addressed, and manipulated, and a single element is 
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already optically dense. Furthermore, the ability to freely convert between phonons and 
photons enables new possibilities for manipulating light through the manipulation of 
sound. Taken together, this raises the possibility that mechanical systems can provide 
a novel, highly configurable on-chip platform for realizing quantum optics and "atomic" 
physics. 
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Appendix A. Equations of motion for optomechanical crystal array 

Here we derive the equations of motion for an array of optomechanical systems coupled 
to a two-way waveguide. Because each element in the array couples independently to 
the waveguide, it suffices here to only consider a single element, from which the result 
for an arbitrary number of elements is easily generalized. 

We model the interaction between the active cavity mode 1 and the waveguide with 
the following Hamiltonian, 

^cav-wg = dk h{ck - uJi)d\^^dR^k dk h{ck + a;i)a^aL,/c 

J— 00 J — 00 

hgV27r f dz6{z - Zj) (a[ (dn^z) + di^z)) + h.c}j . (A.l) 

J— 00 ^ ^ 

Here du^ki^L^k are annihilation operators for left- and right-going waveguide modes of 
wavevector /c, and uoi is the frequency of the cavity mode. For convenience we have 
defined all optical energies relative to cji, and assumed that the waveguide has a linear 
dispersion relation uj{k) = c\k\. The last term on the right describes a point-like coupling 
between the cavity (at position Zj) and the left- and right-going waveguide modes, 
with a strength g. The operator dn^z) physically describes the annihilation of a right- 
going photon at position z and is related to the wavevector annihilation operators by 
dR{z) = ^^00 dke^^^dji^k (with a similar definition for diiz)). Eq. (A.l) resembles 
a standard Hamiltonian used to formulate quantum cavity input-output relations [36], 
properly generalized to the case when the cavity accepts an input from either direction. 
Note that we make the approximation that the left- and right-going waves can be treated 
as separate quantum fields, with modes in each direction running from — oo < A; < oo. 
This allows both the left- and right-going fields to separately satisfy canonical field 
commutation relations, [^^^(z), ai^(z')] = [dL{z)^dL{z^)] = S{z — z')^ while commuting 
with each other. Thus each field contains some unphysical modes {e.g.^ wavevector 
components A; < for the right-going field), but the approximation remains valid as 
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long as there is no process in the system evolution that allows for the population of such 
modes. 

From the Hamiltonian above, one finds the following Heisenberg equation of motion 
for the right-going field, 

V cdi ^ ~dzl ^^^^^ ^ — c — ~ ^^'^^^ ^ ^koaR, (A.2) 

where ko = uji/c. A similar equation holds for a^. The coupling of the cavity mode 
to a continuum of waveguide modes leads to irreversible decay of the cavity at a rate 
/^ex- Below, we will show that is related to the parameters in the Hamiltonian by 
/^ex = ^TTg'^/c. With this identification, one recovers Eq. (5) in the main text. 
The Heisenberg equation of motion for the cavity mode is given by 

= igV27r {aR{zj) + aL{zj)) . (A.3) 

To cast this equation into a more useful form, we first integrate the field equation (A.2) 
across the discontinuity at Zj^ 



o^r{z]~) = aR{zj) + — ^ — ai, (A.4) 



aL{zj ) = aL{zJ) + — - — ai. (A.5) 

We can define aL,in{zj) = aL{zj~) and aR^in{zj) = clr^zJ) as the input fields to the cavity. 
It then follows that 

— ai = z^V27r {an^in + aL,m) (A.6) 

and thus we indeed see that the waveguide induces a cavity decay rate /^ex/2 = I'Kg^ jc. In 
the case where the cavity has an additional intrinsic decay rate /^in, a similar derivation 
holds to connect the intrinsic decay with some corresponding noise input field a^^. From 
these considerations, and including the opto-mechanical coupling, one arrives at Eq. (3) 
in the main text. 



Finally, we consider the equation of motion for the mechanical mode given by 
Eq. (4) in the main text, 

f = -^^ + ^O^al + F^(^). (A.8) 

The zero-mean noise operator F^y must accompany the decay term in the mechanical 
evolution in order to preserve canonical commutation relations of 6, at all times. In 
the case where the decay causes the mechanical motion to return to thermal equilibrium 
with some reservoir at temperature T5, the noise operator has a two-time correlation 
function given by <F^(t)Fj^(f )> = lm{n + \)8{t - f) [37], where n = {e^rn/ksn _ i)-^ 
is the Bose occupation number at the mechanical frequency uOm- 
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Appendix B. Transfer matrix analysis of propagation 

First we derive the reflection and transmission coeflicients for a single element in the case 
of constant opto-mechanical driving amplitude Qju- Given the linearity of the system, it 
suffices to treat Eqs. (A. 2), (A. 7), and (A. 8) as classical equations for this purpose, and 
furthermore to set the noise terms F/v = 0, a^^ = 0. For concreteness, we will consider 
the case of an incident right-going cw field in the waveguide. Upon interaction with the 
opto-mechanical system at Zj^ the total right-going field can be written in the form 

aniz) = e"'=* (e(-^ + zj) + t(5k)e(z - zj)) , (B.l) 

while the left-going field is given by Gl^z) = e~'^^^~'^^^^r{Sk)0{—z + Zj). Here B(z) is 
the unit step function, 5k = ck — uji is the detuning of the input field from the cavity 
resonance, and r, t are the refiection and transmission coefficients for the system. At the 
same time, we look for solutions of the cavity field and mechanical mode of the form 
ai = Ae~^^^^ and h = Be~^^^^ . The coefficients r, t, A, B can be obtained by substituting 
this ansatz into Eqs. (A. 2), (A. 7), and (A. 8). This yields the reflection coeflicient given 
by Eq. (6) in the main text, while the transmission coeflicient is related hy t = 1 + r. 

To calculate propagation through an array of elements, it is convenient to 
introduce a transfer matrix formalism. Speciflcally, the flelds immediately to the right 
of the opto-mechanical element (at z = z'^) can be related to those immediately to the 
left (at z = zj) in terms of a transfer matrix Mom-> 





= Mom ""r'-i - (B-2) 



where 



Mom = -{ , . (B.3) 



On the other hand, free propagation in the waveguide is characterized by the matrix 
Mf, 



aji{z + d) 
ql^z + d) 




(B.4) 



where 



pikd A 



The transfer matrix for an entire system can then be obtained by successively 
multiplying the transfer matrices for a single element and for free propagation together. 
In particular, the transfer matrix for a single "block", defined as interaction with a 
single opto-mechanical element followed by free propagation over a distance d to the 
next opto-mechanical element, is given by Muock = MfMom^ and the propagation over 
blocks is simply characterized by = M^^^j^. 
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Before studying the propagation through the entire array, we first focus on the 
propagation past two blocks, M2 = M^i^^j^. Because we want our device to be highly 
transmitting when the optomechanical coupling is turned on, we choose the spacing d 
between consecutive blocks to be such that k^d = |(2n + 1), where n is an integer. 
Physically, this spaces consecutive elements by an odd multiple of A/4, where A is 
the resonant wavelength, such that the reflections from consecutive elements tend to 
destructively interfere with each other. This can be conflrmed by examining the resulting 
reflection coefficient for the two-block system, 

where M2(i, j) denotes matrix elements of M2. Note that the reflection coefficient is now 
suppressed as a quadratic in the detuning, whereas for a single element r^iHie^5k/2ft^ 
is linear. In the above equation, we have made the simplifying approximation that 
kd^kod^ since in realistic systems the dispersion from free propagation will be negligible 
compared to that arising from interaction with an opto-mechanical element. 

Now we can consider transmission past N/2 pairs of two blocks (i.e., elements 
in total). Because the reflection r2 is quadratic in the detuning, its effect on the total 
transmission is only of order 5^ (because the lowest order contribution is an event where 
the fleld is reflected twice before passing through the system). Thus, up to 0(5f), the 
total transmission coefficient tj^ is just given by ^A^^t^^^, where ^2 = 1/^2(2,2) is the 
transmission coefficient for a two-block system. It is convenient to write tjsf = e^^^^^^ in 
terms of an effective wavevector fceff, which leads to Eq. (7) in the main text. Performing 
a similar analysis for the case where kod = riTT, where reflections from consecutive 
elements interfere constructively, one finds that the bandwidth-delay product for the 
system AcL;rdeiay^l does not improve with the system size. 



Appendix C. Optical noise power 

In this section we derive the optical cooling equation given by Eq. (9) in the main 
text. We begin by considering the system Hamiltonian for a single element (Eq. (1) 
of main text), in the case where the tuning mode 2 is driven on resonance and can be 
approximated by a classical field, a2^a2e~'^^^^ 

Hom = hujid\ai + humb^b + hh(b + fet) (atof2e-^^2' + a^e^^^'^i) • (CI) 

Defining a detuning 5l = uj2 — uji indicating the frequency difference between the two 
cavity modes, we can re- write Eq. (C.l) in a rotating frame, 

Horn = -h5Ld\di + hujjnb^b + hVtm (b + 6"'"^ (cLi + , (C.2) 

where = ha2 (we have re-defined the phases such that is real). In the weak 
driving limit (f^^<A^), the cavity dynamics can be formally eliminated to arrive at 
effective optically-induced cooling equations for the mechanical motion [14, 15]. In 
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particular, the opto-mechanical coupling terms d\b + h.c. and a[b^ + h.c. induce anti- 
Stokes and Stokes scattering, respectively. These processes yield respective optically- 
induced cooling (r_) and heating (r+) rates 

In the case where 5l = — ^m, the cooling process is resonantly enhanced by the cavity, 
yielding a cooling rate Fopt = T_(5l = —uOm) = ^^m/^ given in the main text. Also 

2 

in this case, the optical heating rate is given by r+ = Fopt ^2 ^^5^2 • This leads to the 
net cooling dynamics given by Eq. (9) in the main text, 

dE 

= -Im {Em - f^mflth)-^optEm + ^opt ^2 ^ ^^^^ (Em + fu^m) .(C.4) 

Because the optical cooling process removes phonons from the mechanical system 
via optical photons that leak out of the cavity, one can identify (cui /cum)^ opt Em 
as the amount of optical power that is being leaked by the cavity in the anti- 
Stokes sideband during the cooling process. Similarly, the cavity leaks an amount of 
power {cjJi/(jJm)^opt 2 2 {Em + hUm) in the Stokes sideband. We have ignored this 
contribution in Eq. (10) in the main text, because its large frequency separation {2(jJm) 
from the signal allows it to be filtered out, but otherwise it approximately contributes 
an extra factor of 2 to the last term in Eq. (10). Finally, we remark that the expression 
for Pnoise givcu by Eq. (10) represents an upper bound in that it does not account for the 
possibility that the output spectrum from a single element may exceed the transparency 
bandwidth, which could cause some light to be absorbed within the system after multiple 
reflections and not make it to the end of the waveguide. 

Appendix D. Band structure analysis 
Appendix D.l. Derivation of dispersion relation 

For simplicity, here we work only with the classical equations so that the intrinsic noise 
terms in the Heisenberg-Langevin equations can be ignored. We begin by transforming 
Eqs. (A. 7) and (A. 8) to the Fourier domain. 



0= (iSk-^^b + iQma. (D.2) 

To simplify the notation, we deflne operators at the boundaries of the unit 
cells (immediately to the left of an optomechanical element) given by Cj = —^^/caJl{zy)^ 
dj = —^^/caL{zJ). It is also convenient to re- write the transfer matrix Mom in the form 
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with the parameter l3{5k) given by 

The transfer matrix M^ock describing propagation to the next unit cell can 
subsequently be diagonalized, M^ock = SDS~^^ with the diagonal matrix D given by 

( iKd Q \ 

Physically, this diagonalization corresponds to finding the Bloch wavevectors K{5k) of 
the periodic system. The dispersion relation for the system can be readily obtained 
through the equation 

cos{K{5k)d) = cos(fcrf) - ip{5k) sm{kd). (D.6) 

Writing k in terms of 5k ^ we arrive at kd = uid/c -\- Skd/c. As described previously, 
the desirable operation regime of the system is such that the phase imparted in free 
propagation should be ujid/c = (2n + l)7r/2. For concreteness, we set here ujid/c = 7r/2, 
satisfying this condition. For the frequencies Sk of interest, which easily satisfy the 
condition \6k\ « d/c and ignoring the intrinsic loss /^in, the simple approximate dispersion 
formula 

^o.m = -^(^ (D.7) 

can be found. This dispersion relation yields two bandgaps, which extend from ±HiQ^/2 
and +2f]^//^, in the weakly coupled EIT regime (flrn ^ ^ex)- We therefore have 
three branches in the band structure, with the narrow central branch having a width 
of 4f]^//^ex- This branch has an optically tunable width and yields the slow-light 
propagation. 

The dispersive and lossy properties of the array can also be found by analyzing 
Eq. (D.6) perturbatively. Expanding Eq. (D.6) as a power series in 6k^ we find 

T^/r \ 7 , l^eyi^k ^^ex^in^^^ (2^ex 3/^ex^hi + 12/^ex^m)^fc , ^/r4\ (y\ Q\ 

which agrees with Eq. (7) in the main text. 
Appendix D.2. Fractional Occupation Calculation 

In our system, the Bloch functions are hybrid waves arising from the mixing of optical 
waveguide, optical cavity and mechanical cavity excitations. It is therefore of interest 
to calculate the hybrid or polaritonic properties of these waves, by studying the energy 
distribution of each Bloch mode. 

The number of photons nwG in the waveguide can be found by taking the sum of 
the left- and right-moving photons in a section of the device. Over one unit cell, one 
obtains 
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The relation between this value and the amplitude of the hybrid Bloch wave may be 
found by considering the symmetry transformation used to diagonalize the unit-cell 
transmission matrix. Defining Cj to be the amplitude of the Bloch mode of interest, 
one finds Cj = SuCj and dj = S2iCj^ while from the properties of the symmetry matrix 
'^^ l^jP ^ kjP + l^jP- From here we can deduce the number of excitations in the 
waveguide nwG? the optical cavity and the mechanical cavity for a given Bloch 
wave amplitude: 

nwG = ^\Cj\', (D.IO) 
no = lap (D.ll) 



^ex 

|2 



\cj+d,\' (D.12) 
\sn + S2i\'\Cj\', (D.13) 



rim =\b\' (D.14) 
= J^HttH'- (D.15) 

We then define the fractional occupation in the mechanical mode by nm/{n^^Q + 
no + rim) (with analogous definitions for the other components). These relations were 
used to plot the fractional occupation and colored band diagrams shown in the main 
text. 



Appendix D.3. Pump Input Power 

The technicalities associated with pumping the optomechanical crystal array system are 
subtly distinct from those in its atomic system analogue. Due to the periodic nature of 
the structure and its strongly coupled property (/^ex » ^in), a bandgap in the waveguide 
will arise at the frequency of the "tuning" or pump cavities. This prevents these cavities 
from being resonantly pumped via light propagating in the waveguide. This problem 
may be circumvented by making the pump frequency off-resonant from uj2 (but still 
at uji — uJm)^ and by for example changing the splitting uji — UJ2 to be less than the 
mechanical frequency. Interestingly, the periodic nature of the system also allows one 
to suppress attenuation of the pump beam through the waveguide (which would cause 
inhomogeneous driving of different elements along the array). Here we calculate the 
waveguide input powers required to drive the pump cavities. We then find the effect 
of the cavity dissipation rate tvi^ on the beam intensity, to provide estimates for the 
power drop-off as a function of distance propagated in the system. Finally we note that 
there is a trade-off between required pump intensity and pump power drop-off. In other 
words, by tuning the pump beam to a frequency closer to the pump cavity resonance, 
the required input power is reduced, but the attenuation per unit cell exp(—lm{K}d) is 
increased. 

To calculate the waveguide input power, we start by considering the net photon 
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flux in the right moving direction, 

^R=\cj\'-\dj\\ (D.16) 

Using the properties of the Bloch transformation matrix S and equation (D.13), this 
expression may written in terms of the number of photons in the optical cavity of interest 

where 

The required input power Pin = HuJo^r for the system parameters studied in the main 
text is shown in figure 4. 

To find the attenuation per unit cell exp(— a), with a = lm{K(5k)}d^ of this pump 
beam, we use a perturbative approach similar to that used to find the polaritonic band 
dispersion. By expanding the Bloch- vector as K{5k) = k^^^ + k^~^y5k + . . . and 

using equation (D.6), we find k^~^^ = /^ex/2rf, and k^~^^ = iK^^Ki^/M implying that 
/4:5l. This approximate expression is shown along with the exact calculated 
values for the attenuation in figure 4. 



Appendix E. Implementation in an Optomechanical Crystal 

For the theoretical demonstration of our slow-light scheme, we confine ourselves to a 
simplified model of an optomechanical crystal where only the two-dimensional Maxwell 
equations for TE waves and the equation of elasticity for in-plane deformations of a 
thick slab are taken into account. These equations approximate fairly accurately the 
qualitative characteristics of in-plane optical and mechanical waves in thin slabs, and 
become exact as the slab thickness is increased. In this way, many of the intricacies of 
the design of high-Q photonic crystal cavities (which are treated elsewhere [34]) may 
be ignored, and the basic design principles can be demonstrated in a slightly simplified 
system. 

The two-dimensional optomechanical crystal (2D0MC) system used here utilizes 
the "snowfiake" design [34], which provides large simultaneous photonic and phononic 
bandgaps in frequency. Here we choose to use optical wavelengths in the telecom band, 
z.e., corresponding to a free-space wavelength of A ^ 1.5 /im. For this wavelength, we 
found that the crystal characterized by a lattice constant a = 400 nm, snowflake radius 
r = 168 nm, and width w = 60 nm, shown in figure 2, should work well. 



Appendix E.l. Optical and Mechanical Cavities 

Appendix E.1.1. Single Cavity System We begin our design by focusing on the creation 
of a single optomechanical cavity on the 2D0MC, with one relevant optical and 
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mechanical mode. This cavity is formed by creating a point defect, consisting of 
two adjacent removed holes (a so-called "L2" cavity). We calculate the optical and 
mechanical spectra of this cavity using COMSOL, a commercial FEM package, and find 
a discrete set of confined modes. Of these, one optical and one mechanical mode were 
chosen, exhibiting the most promising value of the opto-mechanical coupling strength 
g (see below for calculation). These modes are shown in figure 2(b) and (c), and were 
found to have frequencies = 11.2 GHz and Uo = 199 THz, respectively. 

Appendix E.l. 2. Double Cavity System From here we move to designing the nearly- 
degenerate double optical cavity system with large cross-coupling rates. As two separate 
L2 cavities are brought close to one another, their interaction causes the formation of 
even and odd optical and mechanical super-modes with splittings in the optical and 
mechanical frequencies. This splitting may be tuned by changing the spacing between 
the cavities. We take the even and odd optical modes of this two-cavity system as our 
optical resonances at (j02 and uoi. 

Appendix E.l. 3. Optomechanical Coupling Rates The optomechanical coupling arises 
from a shift in the optical frequency caused by a mechanical deformation. Our 
Hamiltonian for the single cavity system can then be written as 

H = tyjj{x)a^a + humb^b, (E.l) 

where x = Xzpf(&^ + b) is the quantized displacement of the mechanical mode, and Xzff 
is the characteristic per-phonon displacement amplitude. The deformation-dependent 
frequency uj{x) may be calculated to first order in x using a variant of the Feynman- 
Hellman perturbation theory, the Johnson perturbation theory [38], which has been 
used successfully in the past to model optomechanical crystal cavities [9, 39]. The 
Hamiltonian is then given to first order by 

H = huJoCL^a + huomb^b + hg{b^ + b)a^a, (E.2) 

where uoq is the optical mode frequency in absence of deformation and 



-^^d/(Q.n) (A6|Ellp-A(6-^)|D^p) 



9 = ; ^ (E.3) 



Here E, D and Q are the optical mode electric field, optical mode displacement field 
and mechanical mode displacement field, respectively, ds is the thickness of the slab, 
and e(r) is the dielectric constant. 

These concepts can be extended to optically multi-mode systems, represented by 
the Hamiltonian 

H = hj^ujo^ialdi + fihJmb'^b + ^J^gijib"^ + b)a|a^-, (E.4) 
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where now the cross-couphng rates can be calculated by the following expression: 
^ . r^^dl (Q . n) (AeEf* . Ell - Aie-^)Bi* • D^^) 



2uJmds 



^dAp\Q\^^dAe\Ei\^^dAe\Ej 



We denote this expression for convenience as gij = (E^|Q|Ej). 

For the modes of the L2 cavity shown in figures 2(b) and (c), the optomechanical 
coupling was calculated to be (E|Q|E)/27r = 489 kHz for silicon. When two cavities 
are brought in the vicinity of each other, super-modes form as shown in figure 2(a). We 
denote the symmetric (+) and antisymmetric (— ) combinations by E+ and Q+. These 
modes can be written in terms of the modes localized at cavity 1 and 2, 

El + Eo Qi + Q9 .X 

E± = and = ^^^^^ (E.6) 

By symmetry, the only non- vanishing coupling term involving the Q_ 
(antisymmetric mechanical) mode is <E+|Q_|E_). Assuming that the two cavities are 
sufficiently separated, we can approximate (E+|Q_|E_) ^ (E|Q|E)/V2. For the super- 
modes of interest, <E+|Q_|E_>/27r = h/27r = 346 kHz. 



Appendix E. 2. Properties of Snowflake Crystal Waveguides 

A line defect on an optomechanical crystal acts as a waveguide for light [40, 41]. Here, 
the line defects used consist of a removed row of holes, with the rows above and below 
shifted towards one another by a distance such that the distance between the centers 
of the snowfiakes across the line defect is ^/3a — 2W (see figure El (a)). The waveguide 
was designed such that mechanically, it would have no bands resonant with the cavity 
frequency (see figure El(b)) and would therefore have no effect on the mechanical Q 
factors. Optically, it was designed have a single band crossing the cavity frequency (see 
figure El(c)) and would therefore serve as the single-mode optical waveguide required 
by the proposal. The band structure of the mechanical waveguide was calculated using 
COMSOL [42], while for the optical simulations, MPB [43] was used. 



Appendix E. 3. Cavity- Waveguide Coupling 

By bringing the optical waveguide near our cavity, the guided modes of the line-defect 
are evanescently coupled to the cavity mode, and a coupling between the two may be 
induced, as shown in figure E2(a). Control over this coupling rate is achieved at a coarse 
level by changing the distance between the cavity and waveguide, i.e., the number of 
unit cells between them. We found a distance of 6 rows to be sufficient in placing our 
coupling rate A^ex in a desirable range. At this point, a fine tuning of the coupling rate 
may be accomplished by adjustment of the waveguide width parameter W^ described 
previously. The achievable values of 

Avgx are plotted against W in figure E2(b). For the 

final design, W = 0.135a is used. 
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Figure El. (a) The snowfiake waveguide consists of the usual snowfiake lattice, with 
a removed row, and adjusted waveguide width. The guiding region is shaded red for 
clarity, (b) The mechanical band structure exhibits a full phononic bandgap at the 
resonance frequency of the mechanical mode (i^rn)- (c) The optical band structure 
exhibits a single band that passes through the resonance frequency uq of the optical 
mode. The waveguide acts as a single- mode optical waveguide with field pattern 
shown in (d). The Hz component of the guided optical mode of opposite symmetry is 
shown in (e) for completeness. 



To simulate this coupling rate, we performed finite-element simulations using 
COMSOL where we placed the waveguide near our cavity, and placed absorbing 
boundaries at the ends of the waveguide away from the cavity. The resulting time- 
averaged Poynting vector (S)^ = |E x H*|/2 is plotted in figure E2(a), showing how the 
power flows out of the system. 

Appendix E.Jf.- Estimates for Thermoelastic Damping 

The achievable storage times of our system are determined by the lifetimes of the 
mechanical resonances. Since we use a phononic crystal, all clamping losses have been 
eliminated. However, other fundamental sources of mechanical dissipation remain, and 
here we provide estimates for one of these, the component due to thermoelastic damping 
(TED) [44, 45]. 

Using the COMSOL finite-element solver [42], we solved the coupled thermal and 
mechanical equations for this system [46] . In these simulations the change in the thermal 
conductivity and heat capacity of silicon with temperature were taken into account. The 
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Figure E2. (a) A plot of the magnitude of the time- averaged Poynting vector, 
|(S)^| (in arbitrary units), shows the leakage of photons out of the double-cavity 
system. This induces a loss rate on the optical modes, which was used to calculate the 
extrinsic coupling rate A^ex, plotted in (b) for various values of W/a. (c) Plot of the 
mechanical quality factor (5m,TED due to thermoelastic damping, as a function of the 
ambient temperature Tq. (d) The time-harmonic component of the temperature field 
AT(r) = T — To is plotted at various times during a mechanical oscillation period. 



TED-limited quality factors, (5m,TED are plotted in figure E2(d). In these simulations, 
we see that for the mode simulated, (5m,TED surpasses 10^ at bath temperatures of 
To < 5 K. To illustrate some representative results of these simulations, we have plotted 
the change in temperature field AT(r) from the ambient temperature versus the phase 
of the mechanical oscillation in figure E2(d). At cj^t = 7r/2, there are variations in 
temperature despite the displacement field Q being uniformly at this time. This 
shows that at these frequencies, the temperature does not follow adiabatically the 
displacement. 

Appendix E. 5. Estimate for Optical Pump Heating 

As mentioned in the main text, in a realistic setting the optomechanical driving 
amplitude itself will be coupled to the bath temperature through absorption of 
optical pump photons in the tuning cavities. This optical pump heating of the structure 
is important in estimating the practical limits of the optomechanical system for quantum 
applications where thermal noise impacts system performance. As a realistic model for 



Slowing and stopping light using an optomechanical crystal array 



26 



the bath temperature in our proposed Silicon optomechanical crystal array, we take 
= To + x^^i where Tq is the base temperature and x is a temperature coefficient that 
describes the temperature rise in each cavity per stored cavity photon due to optical 
absorption. Our estimate of x foi" ^ thin-film Silicon photonic crystal structure is as 
follows. The absorbed power for |ap photons stored in a cavity is given simply by 
Pioss = huoi^il^^f^ ^ where Uo is the resonance frequency and HCi the optical (intrinsic) 
linewidth of the cavity. If we assume all of this power is being converted to heat, the 
change of temperature is AT = PiossRth, where Rth is the effective thermal resistance 
of the Silicon structure. There are a number of sources in the literature for Rth in 
relevant photonic crystal geometries [47, 48]. We choose here to use the value for a 
two-dimensional crystal system in Silicon, Rth ^ 2.7 x 10^ K/W, which yields a per 
photon temperature rise of x^2 /iK assuming an intrinsic loss rate of ^ 4 x 10^ rad/s 
(Q, ^ 3 X 10^). 
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